Low density expansion and isospin dependence of nuclear energy functional: 
comparison between relativistic and Skyrme models 
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In the present work we take the non relativistic limit of relativistic models and compare the 
obtained functionals with the usual Skyrme parametrization. Relativistic models with both constant 
couplings and with density dependent couplings are considered. While some models present very 
good results already at the lowest order in the density, models with non-linear terms only reproduce 
the energy functional if higher order terms are taken into account in the expansion. 
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I. INTRODUCTION 

Both non-relativistic and relativistic phenomenological 
nuclear models are nowadays used to describe success- 
fully nuclear matter and finite nuclei within a density 
functional formalism. One of the most commonly used 
non-relativistic models, the Skyrme force model, has been 
extensively used in the literature since the work of Vau- 
therin and Brink [1]. The Skyrme effective interaction 
has been improved since its first version [2] in order to 
account for the different properties of nuclei, not only 
along the stability line, but also for exotic nuclei from 
the proton to the neutron drip line and to yield a good 
description of neutron stars [3]. Along this line, many 
different versions of the Skyrme effective force have been 
developed and tested and the importance of the asymme- 
try discussed [3, 4]. In [5], the results of the variational 
microscopic calculations, in which many-body and some 
relativistic corrections were included, were parametrized. 
In [4] the Skyrme model parameters were chosen so as to 
fit the EoS given in [5]. In [6] 87 different parametriza- 
tions, which give similar results for finite nuclei exper- 
imental observables, were checked against neutron star 
properties and 60 of them were ruled out. Some of the 
non-relativistic models include also a three-body force [7] 
in order to solve the causality problem, reproduce satu- 
ration properties of nuclear matter and improve the de- 
scription of the symmetry energy. 

Relativistic models, on the other hand, are advanta- 
geous if high density matter, such as the one existing in 
compact stars, is described, and in particular no lack of 
causality arises as density increases. As seen in [6], neu- 
tron star properties depend on the correct choice of the 
equation of state (EoS) and the same problem appears 
within relativistic models, which can be parametrized in 
different ways, all of them giving similar results for finite 
nuclei and nuclear matter or neutron matter. Once the 
EoS is extrapolated to high densities, quantities as sym- 
metry energy, for example, depend a lot on the param- 
eter set chosen [8, 91. At subsaturation densities, there 



is an unstable region, which varies a lot and present dif- 
ferent behaviors according to the parametrization chosen 
[10-12]. Some relativistic models introduce the density 
dependence through the couplings of baryons to mesons 
[13-15] and once again the instability region at low den- 
sities and the EoS at high densities are sensible to the 
model considered [11, 15]. We point out that the satura- 
tion mechanism of relativistic and non-relativistic mod- 
els is different. For the first ones saturation is attained 
due to the relativistic quenching of the scalar field. On 
the other hand, non-relativistic models have to introduce 
three-body repulsive interactions in order to describe sat- 
uration correctly. 

At very low densities both, the relativistic and the non- 
relativistic approaches predict a liquid-gas phase transi- 
tion region for nuclear matter leading, for dense star mat- 
ter to a non-homogeneous phase commonly named pasta 
phase, formed by a competition between the long-range 
Coulomb repulsion and the short-range nuclear attrac- 
tion [16]. 

In the past some attempts have already been made in 
order to compare nuclear matter and finite nuclei prop- 
erties obtained both with relativistic and non-relativistic 
models [11, 17, 18] but there is no clear or obvious expla- 
nations for the differences. In the following we compare 
the Skyrme effective force with relativistic nuclear mean- 
field models at subsaturation densities with the goal to 
directly compare the energy functional. We expect that 
the same physics should be contained in both approaches 
at these densities and therefore it maybe fruitful to com- 
pare both types of models in this range of densities. 

We start with a brief review of the Skyrme 
parametrizations of the nuclear energy density func- 
tional. Next we consider relativistic mean-field models 
with constant couplings, namely the initially proposed 
parametrizations including only linear terms for the me- 
son contributions [19] and parametrizations with non- 
linear terms [20-22]. Then we extend our investigation 
to density dependent models [10, 13, 14, 23, 24]. We 
compare the non-relativistic limit of their binding ener- 
gies with the binding energy functional obtained with 
the non-relativistic Skyrme model. Various levels of the 
approximation are then discussed. 
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II. SKYRME FUNCTIONAL 



The non relativistic Skyrme energy functional is de- 
fined by 



Bsky 



K, + Ho + H 3 + Ti-eS, 



where K is the kinetic-energy density, Ho a density- 
independent two-body term, H 3 a density-dependent 
term, and H c s a momentum-dependent term: 



K 

Ho 
H 3 

H e f[ 



2M' 

C oP 2 + D oP l 
C 3 p a+2 + D 3 p° pi 
C e ffpr + D cS p 3 r 3 , 



(1) 

(2) 
(3) 
(4) 



where p = p p + p n is the total baryonic density, p p and p n 
being the proton and the neutron densities respectively, 
A3 = Pp-pn, the isovector density, r = J2i= P . n T ^ tnc to " 
tal kinetic density, T{ — k\. /5tt 2 being the kinetic density 
of each type of particles i with a Fermi momentum 
and t 3 = t p — r n the isovector kinetic term. We intro- 
duce the following quantities: B\ the symmetric matter 
potential energy, B 3 the potential part of the symmetry 
energy and M* the proton (i = p) or neutron (i = n) 
effective masses, 



Bi{p) 
B 3 (p) 
Mf 1 

such that 

-^Skyrme 



Cop 2 + C 3P ° +2 , 
Do + D 3 p°, 

M" 1 + 2C eS p + t 3i 2D cS p 3 , 



= E 



Tj 

2M; 



■+B 1 (p)+p 2 3 B 3 (p). 



(5) 
(6) 
(7) 



(8) 



The coefficients Ci and Di, associated respectively 
with the symmetry and asymmetry contributions, are lin- 
ear combinations of the traditional Skyrme parameters : 

Co = 3to/8 

Do = -t (2x + l)/8 

C 3 = t 3 /16 

D 3 = -t 3 (2x 3 + l)/48 

C cS = [3t 1 +t 2 (4x 2 + 5)}/W 

D c g = [t 2 (2x 2 + l)-ti(2xi + l)]/16 

and a parametrizes the density dependent term. The a 
exponent has a direct effect on the incompressibility. A 
decrease in a generates a lower value for the incompress- 
ibility [3]. Using pi = kp./3-K 2 the kinetic terms can be 
directly related to the respective density 

n - a(p,) 5/3 

where we have introduced a = 3 5 / 3 tt 4 / 3 /5. Looking at 
small asymmetries p p — (p ± p 3 ) /2 we get 



t 3 = 



> 1 p 5/3 + Ip-^pI 

5a 



P 2/3 P3 



(9) 
(10) 



Then it is easy to decompose the kinetic energy density 
ICdf =K + HeS into a symmetric matter component and 
a symmetry kinetic energy 



leading to 

^cffi 



K-ett = ^cffi + P 3 fc e ff 3 



ICcS 3 — a-p 



'2M 

-1/3 



2M 



+ (C eff + 3D eS )p 



The isoscalar part of the effective mass C c s directly con- 
tributes to /Ceffi but also contributes to the symmetry 
energy because of the asymmetry in the Fermi energy. In- 
deed as in the Fermi gas model the two first terms of /C e ff 3 
can be recasted as IC e s 3 — 5/0^ /9,o 2 . The last term is 
an additional term coming from the isovector part of the 
effective mass. These various contributions to the kinetic 
energy play an important role both in the binding energy 
and saturation properties of symmetric matter and in the 
isospin dependence usually discussed in terms of symme- 
try energy. When comparing with relativistic approaches 
the discussion of the kinetic part becomes even more im- 
portant since the role played by the effective mass and 
interaction are known to be in general rather different. In 
Table I we present the nuclear matter saturation proper- 
ties obtained with the Skyrme forces used in the present 
work. We consider a conventional Skyrme interaction 
SIII, one of the recent Skyrme-Lyon interactions designed 
to describe neutron-rich matter SLy230a [3] , the NRAPR 
parametrization which stands for the Skyrme interaction 
parameters obtained from a fitting to the EoS of a micro- 
scopic model [4, 5] and the LNS parametrization which 
refers to a recently Skyrme- like parametrization proposed 
in the framework of the Brucckner-Hartree-Fock approx- 
imation for nuclear matter [25]. 



III. RELATIVISTIC APPROACHES 

A. The lagrangian 

In this section we consider four mean-field relativis- 
tic models, which we denote by QHD-II [19], NL3 [20], 
TM1 [21] and NLS [22], with constant coupling param- 
eters described by the Lagrangian density of the linear 
and non-linear Walecka models (NLWM), with the pos- 
sible inclusion of the S mesons, given by: 



"NLWM 



= A +C a +C ul +C p +C 5 , (11) 
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where the nucleon Lagrangian reads 



with 



iD" = id" -g v V» -^T-b», 
M* = M - g s <j>- g s r -6. 



(12) 

(13) 
(14) 



The isoscalar part is associated with the scalar sigma 
(a) field, <f), and the vector omega (u) field, V^, while 
the isospin dependence comes from the isovector-scalar 
delta (S) field, 5 l , and the isovector- vector rho (p) field, 
b p (where p is the 4 dimensional space-time indices and 
i the 3D isospin direction indices). The associated La- 
grangians are 



+ \ (d^O^ - m 2 ^ 2 ) - I K 3 - iA0 4 , 
+ l -{dJd»S-m 2 s 5 2 ), 



Ha = 

r - 

Cs = 
H P = 

where M „ = d M V„ - , B^ = dj> v - d„6 M - g^b^ x 
&„), and gj and are respectively the coupling con- 
stants of the mesons j = s, v, p, S with the nucleons and 
their masses. Self-interacting terms for the cr-meson are 
also included in the three parametrizations, k and A de- 
noting the corresponding coupling constants. The w- 
meson self-interacting term, with the £ coupling constant, 
is present in the TM1 parametrization. In the above La- 
grangian density f is the isospin operator. The terms in- 
volving the S meson are only present in the NL-S model 
and non-linear terms are not present in the simplest ver- 
sion of the model QHD-II. In Table II we present the 
nuclear matter saturation properties obtained with the 
relativistic models used in the present work. 



B. Equilibrium 

At equilibrium the fermion distribution is a Fcrmi- 
Dirac distribution function for nucleons (+) an antin- 
ucleons (-) 



/ t± = l/{l + eX p((e*(p)T^)/T)} 



(15) 

where e* = \J\> 2 + M* 2 with the effective mass 

M*=M- g s cb Q - r 34 gs 5 3 , (16) 

In the equilibrium, the effective chemical potentials are 
defined by 



n = Hi - g v Vo 



2 



T 3 i b , 



(17) 



where T 3 i = ±1 is the isospin projection for the protons 
and neutrons respectively. At zero temperature the dis- 
tribution function reduces to a simple step function 



fi = 0(pl-p 2 ), 



with 



M* 



(18) 



(19) 



Let us introduce <po, Vo, S 3 and bo the values of the 
scalar (<t), the vector (u>), isovector scalar (6) and the 
isovector vector (p) fields, obtained from the meson equa- 
tions of motion, considered as static and uniform classical 
fields 



m 2 s (j)o + ^ 




= g s p s , 


(20) 


m 2 y 


+ JfffX 


= g v p, 


(21) 




m 2 6 3 


= gsPs3, 


(22) 




m 2 p b 


= yP3, 


(23) 



where the baryonic and scalar densities read: 

d 3 p 



p 4 = 2 



d 3 p 

(2^) 3 



(2tt) : 



M* 



*2 



-ApI -p 2 ), 



(24) 



(25) 



and p = p p +p n , p 3 = Pp-pn , Ps = Psp+Psn , Ps3 = Psp~ 

p sn . The above ensemble of equations for the effective 
mass (16), the fields (20)-(23) and the densities (24)-(25) 
related to the Fermi momentum (19) defines entirely the 
state in a self consistent way. In the sequel we introduce 
the constants c p = g 2 /8m 2 and c a = </ 2 /2m 2 , with a = 
s, v, S. 

The energy density is then given : 



E = K P + K n + E ir + E u + Es + E p , 



with 



Ki 

E a = 

E u = 

E 5 = 

E n = 



1 



1 



2 

^4 + ^4 + ^\4, 
^v 2 + 

m S c-2 2 

— S 3 = c s p s3 , 
c pP 3 , 



(26) 

(27) 

(28) 

(29) 

(30) 
(31) 



Having solved the field's self-consistent equations entirely 
defines the system energy functional. 
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In order to compare the relativistic approaches with 
the non relativistic one we also look at the binding energy 
density defined as 

B = E — Mp. (32) 
Isolating the kinetic contribution 

Ti = Ki- M*p, (33) 
we can define an interaction part 

B exact (p,p 3 ) = B-J2 T h (34) 

i 

which can be decomposed into an isoscalar part 

-Bi oxact (p) = B cxact {p, p 3 = 0) (35) 
and isovector potential energy contributions 

D , n B exact (p, p 3 ) — -Blexact(p) ,„„n 

i>3cxact(P, P3) = 2 ' V-JOj 

P3 

Using these notations the binding energy from relativistic 
approaches can be recasted as in the Skyrme case 

B = ^TT t + Blexact(p) + B 3cxact (p, p 3 ) p\. (37) 

i 

We also separate the effective mass into the isoscalar 
and the isovector channels and write 




0.05 0.1 0.15 0.2 0.25 0.3 

p (fm -3) 



FIG. 1: Symmetric matter: a) Binding energy; b) Isoscalar 
contribution of the interaction divided by p 2 , Bicxact/p 2 ; c) 
Effective mass for several relativistic and non-relativistic mod- 
els. 



M* =M + Mi + T 3l p 3 M 3 . 



(38) 



gives 



IV. DIRECT COMPARISON OF ENERGY 
FUNCTIONAL 

The first idea to make a bridge between the classical 
and relativistic models is to directly compare the various 
terms entering in the associated energy functionals. In- 
deed, considering the Kohn-Sham theorem, models lead- 
ing to the same energy functional are strictly equivalent 
as far as physical results are concerned. However, in the 
previous sections, we have seen that for spin saturated 
static matter the Skyrme energy is a functional of the 
isoscalar and isovector particle and kinetic densities. The 
case of the relativistic models is more complex since the 
proton and neutron baryonic and scalar densities appear 
as well as the relativistic kinetic energy density. This 
makes a direct comparison of the functional impossible. 
In this paper we thus adopt two strategies: i) In the 
next section we develop a non-relativistic approximation 
to these densities in order to reduce them to the Skyrme 
like local densities. Then a direct comparison of the func- 
tionals at low densities is possible and will teach us a lot 
about the connections between those two classes of mod- 
els. A low density expansion of the energy per particle 



E/p 



M+ ^ 



3k p 

10M~56M 3 



+ 



9]_ P_ 

2 M 



m 



10M 



.+ (-h- 

+ \2to2 

( 9sP \ 2 
\mlMJ 



9 s \ 



10M 



+0 



' 9 8 p y 



where the terms of the type y ^i M J , n = 1,2,3,... 
have origin on the expansion of the scalar density and 
are equivalent to n-body terms of the Skyrme interac- 
tion generally described by fractionary exponents, ii) 
Secondly, we look at the contributions of the various en- 
ergy terms to the nuclear matter EoS thus reducing the 
functionals of both models to particle densities. With 
this reduction, exact results for the various models can 
be directly compared. 

Let us first focus on the properties of symmetric mat- 
ter. In Figs, la), lb) and lc) we show respectively the 
binding energy, the isoscalar contribution of the interac- 
tion -Bioxact and the effective mass for several relativistic 
and non-relativistic models. We include the results of the 
Walecka (QHD-II) model just for comparison although 
we know it does not reproduce well the properties of nu- 
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clear matter at saturation. In particular, it has a very 
high compressibility. A comparison of the binding ener- 
gies plotted in Fig. la) shows that i) although the satu- 
ration point is quite close in energy and densities for all 
models, there is some discrepancy which is related to the 
way the paramctrizations were fitted; ii) relativistic and 
non-relativistic models show different behaviors above 
an below saturation density: at subsaturation densities 
relativistic models show generally more binding (except 
for the Walecka parametrization) and above saturation 
density the binding energy decreases faster with density 
within these models. The parametrization with density 
dependent coefficients (TW) has a behavior closer to the 
other non-relativistic models; iii) the second derivatives 
of the curves are associated with the incompressibility 
and can be very different. This is just a reflex of the 
values tabulated for the compressibilities of the different 
models in Tables I and II. 

We next compare the isoscalar term of the interaction. 
In order to separate the parabolic contribution we repre- 
sent this term divided by p 2 . The Walecka parametriza- 
tion shows an almost independent behavior on the den- 
sity because this model has no non-linear terms on the 
CT-mcson. Non-relativistic and relativistic models show 
similar behaviors although the last ones do not grow so 
fast with density. This tendency will be compensated 
by the density dependence of the effective mass in the 
two types of models. In fact, the effective mass decreases 
faster in the relativistic models as compared with non- 
relativistic models, as seen in Fig. lc). Smaller effective 
masses give rise to larger kinetic energy contributions so 
that the kinetic energy contribution to the total energy is 
more important in relativistic models. We should stress, 
however that the effective masses in relativistic and non- 
relativistic models have different meanings. For the first 
type the effective mass includes the contribution of the 
nucleon scalar self-energy, while for the second type the 
effective mass reflects the momentum dependence of the 
single particle energy. We point out that the isoscalar 
channel of the Skyrme parametrizations SLy230a and 
LNS are very similar. These two parametrizations show 
differences in the isovector channel. 

Let us now turn to the isospin dependence. We focus 
the discussion on the symmetry energy, the exact isovec- 
tor term of the interaction i? 3 [Eq. (36)] and the isospin 
channel of the effective mass M3. We include the rela- 
tivistic models with non-linear terms on the p — oj mesons 
which allow a more flexible description of the isovector 
channel. From Fig. 2a), one observes that the slopes of 
the symmetry energy are somewhat different when con- 
sidering a relativistic or a non-relativistic model, TW 
being the only relativistic model that shows a behaviour 
similar to the non-relativistic ones. The other relativis- 
tic models show an approximately linear dependence of 
the symmetry energy on the density. It has recently been 
discussed in [26] that there is a linear correlation between 
the neutron skin thickness and the slope of the symme- 
try energy for non-relativistic models and this fact gives 




0.05 0.1 0.15 0.2 

P (fm -3) 



FIG. 2: Isospin channel: a) symmetry energy; b) exact B3; c) 
isospin channel of effective mass M3. 



a very strong constraint on the density dependence of the 
nuclear symmetry energy and consequently on the EoS 
as well. In [27] it was found that the same kind of corre- 
lation exists for density dependent relativistic models. 

The isovector channel of the interaction, Fig. 2, also 
shows a large discrepancy between the different models 
and also between models within the same framework: in 
particular the SLy230a shows a behavior very different 
for the other Skyrme parametrizations, with a very fast 
decrease with density becoming even negative for p > 
0.15 fm -3 . The relativistic models with the non-linear 
lu— p terms or with density dependent couplings including 
the 8 meson also become negative at large densities but 
this effect is not so pronnounced. 

In relativistic models there is a proton-neutron mass 
splitting only if the scalar isovector (5-meson is included. 
This occurs for the DDH<5 and NL<5 parametrizations we 
consider. As we see from Fig 2c), M 3 = (M* - M*)/p 3 
is negative for these models which corresponds to a 
M* < M* in neutron rich nuclear matter. A similar 
behavior is predicted by the Skyrme interaction SLy230a 
but an opposite behavior is obtained with the LNS and 
NRAPR parametrizations of the Skyrme interaction. A 
discussion of the isospin dependence of the effective mass 
has been done in [28] . The proton-neutron mass splitting 
is a present topic of discussion and the forecoming exper- 
iments with radioactive beams will allow the clarification 
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of this point. 

From the figures just discussed we conclude that the 
isoscalar channel already shows some discrepancy be- 
tween the different models, mainly between the relativis- 
tic and the non-relativistic ones. However the larger dis- 
crepancies occur for the isovector channel. In this case 
there is not even a general common trend between the 
models of each class. 



V. NON-RELATIVISTIC APPROXIMATION TO 
THE ENERGY FUNCTIONAL 

In what follows we perform a non-relativistic expansion 
of the energy density [Eq. (26)] in order to make a com- 
parison with a classical approximation like the Skyrme 
one. This approximation is based on the expansion 



P 2 + M, 



*2 



Mf- 



2m; 



(39) 



In order to be valid the Fermi momentum should be small 
compared to the nuclcon effective mass meaning that 
such an expansion is restricted to low densities. Then 
we have the scalar density and the relativistic kinetic 
energy written as functional of the particle and kinetic 
densities: 



2M* 



(40) 



In the same spirit, the a field is directly deduced from 
Eq. (20) leading to 



sps/m 2 s 



(44) 



Using the expansion (40) of the scalar density p s , the a 
field can be explicitly written as a functional of p and t 



g s (j> = 2c s (p 



2M* 



(45) 



Since in this model the a field is the only contribution to 
the effective mass, using (16), the effective mass is given 
by: 



M* = M -2c s (p 



2M* 



(46) 



This equation can be solved iteratively replacing 1/M* 2 
on the rhs of equation (46) by M~ 2 + Ac s pM~ 3 . The 
different terms of the energy (26) arc now all functional 
of the particle densities and kinetic densities. 



E n 



M -2c. [p- 



2M* 



2M*' 



c s p 



c v p 



M* 



m 



PU2 



(47) 

(48) 
(49) 
(50) 



Ki ~ M* Pi + 



2M* 



(41) 



Using these two approximations we have reduced the rel- 
ativistic densities to the standard densities appearing in 
the Skyrme approaches. Solving now the field equations 
it becomes easy to recast the relativistic energy func- 
tional in a Skyrme form for the considered case, depend- 
ing only on particles and kinetic densities. Since this 
equivalence is obtained at the first order in p 2 /M 2 we 
obtain analytical expressions for the fields at the same 
level of approximation. 



A. QHD-II 

In order to test the non-relativistic approximation 
of the relativistic models we first consider the QHD-II 
parametrization [19]. In this case the non- linear terms 
are zero (n = 0, A = 0, and £ = 0) and the 5 meson is 
not present (c$ = 0). 

Then, the u field is directly deduced from Eq. (21) 



V = g v p/m 2 v , 
and so is the p field from Eq. (23) 

9p 



bo 



-P3/m p . 



(42) 



(43) 



and the energy functional can be analytically compared 
with the Skyrme one, the differences being both in the 
functional dependence and in the coefficients. Then 
grouping the Ki and E a terms helps to recognize the 
different terms of the Skyrme functional: 



i—p.n 



with 



EL = 



2M*' 

-CsP 2 , 



where the effective mass is approximated by Eq. (46). 

We can now look at the binding energy density given 
by Eq.(32) which becomes 

-B n on-rel. = ^ Ti + E' a + E u + E p . (51) 
i—p.n 

In order to test this Skyrme-like approximation, we 
rewrite the functional B n on-roi. in the form of Eq. (8) 
and obtain 



Bi{p) = ( c v -c s )p 2 
B 3 (p) = c p . 



(52) 
(53) 
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In Figs. 4a) we compare the approximate binding energy, 
-B n on-rci ., Eq. (51), with its exact value, Eq. (32), the 
quantities -Bi(p), eq (52) with -E>i oxact (p) defined in Eq. 
(35) and M*, Eq. (46) with M c * xact = M - g s (j> a . In all 
cases we can see that the non-relativistic approximation 
is a good approximation for low densities when the Fermi 
momentum is not too high. Thus, we can safely compare 
the functional obtained in this limit with the Skyrme one. 

We can see that the above non-relativistic limit of this 
first simple relativistic model (QHD-II) leads to a Skyrme 
functional as far as the density dependence is concerned 
(see Eqs. (5) and (6) ), with a simple two body force (C3 
and D 3 are both zero) . As already referred, we point out 
that the three-body term in the Skyrme parametrization 
is essential to get saturation. For the simple Walecka 
model the quenching of the scalar field will play a similar 
role. The effective mass is similar but not identical since 
the Skyrme parametrization Eq. (7) applies to M* 
while the relativistic models introduce directly M* (see 
Eq. (46)). Moreover the relativistic models in their non- 
relativistic limit present a richer parametrization of the 
effective mass with non linear contributions in the kinetic 
energy density. However, the leading orders are similar. 

It should be noticed that, for this simple model, with- 
out a S field the isospin dependence of the effective mass is 
not introduced. In the next section we study the isospin 
dependence of the mass. 

In order to introduce non-trivial density dependence, 
we also consider models with interacting fields such as 
the non-linearities in the a models or the non-linear up 
and wp couplings (see section VC). An alternative way 
is to directly introduce density dependent coupling pa- 
rameters as discussed in section VD. 



B. S and non-linearities in the a and uj fields 



the proton and neutron effective masses, Eq. (40), which 
is a function of the a field, Eq. (16). In the case of small 
non-linear terms, the a field can be perturbatively solved. 

Indeed, in a simple and crude approximation, the k 
and A terms are assumed to be small, and then from Eq. 
(20) the leading term becomes: 

— T T 

g s <f> = 2c s (p Sp +p Sn ) = 2c s (p - - 2]^)- ( 55 ) 

This approximation can thus be introduced in the non- 
linear terms to get 0o = 4>o + d(j) with, at the lowest order 
in the non linear terms, 

*-5^*-i4* <56) 

This leads to 

c s 

4>o = 2 — p - A CT - Ako, 
9s 

with the non-linear and kinetic contributions 

mjgj 3 mjgj 

Ako ~ Vs { W w l 

The same strategy can be used to solve the w field 
either directly from its relation to the baryonic density 
(cf. Eq. (21)) or assuming that at low density the non- 
linearities are small and thus solving Eq. (21) iteratively 
from the linear solution Vo = g v p/m 2 . Then at the first 
iteration we get 



In this section, we cure the two problems encountered 
in the simple model discussed above introducing the S 
field in order to produce an isospin dependent mass split- 
ting between protons and neutrons and a non-linear a 
and lo field to modify the density dependence of the en- 
ergy functional. 

The first step is to solve the field equations. Let us 
first start with the 8 field which can be easily solved in- 
troducing the approximate scalar field (40) in the field 
equation (22). We obtain for the S 3 field 



with 



2— [Ps 

gs 



A/f3 = 



PsJ = 2 p 3 

gs 



c <5 / 7~p T n . 



*K3, 



(54) 



Turning now to the a field we have to solve the self- 
consistent problem relating 4> to the proton and neutron 
scalar fields (20) which can be computed as a function of 



rat 6 m° 



(57) 



The expressions for </> and £3 can now be used to ob- 
tain the non-relativistic approximation for the effective 
mass in the form of Eq. (38), with 

Mi = -2c s p + g s A a +g s A Ka , (58) 
M 3 = -2c 5 +g s A K: Jp 3 . (59) 

This set of coupled equations can be solved directly or 
iteratively. At the first iteration the (M*)~ 2 reads 

(M*)- 2 = M- 2 + 4c s pM- 3 + 4c S T 3lP3 M- 3 . (60) 

Having solved the equations for the fields and for the 
effective masses in terms of the particle and kinetic den- 
sity we can now study the energy functional. 

In the present case, the different terms of the energy 
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read 



Ki = p z (M - 2c s p - t 3i 2c s p 3 + g s A a 

T' 

+g s A KQ + gsT 3i A K3 ) + 



2M* 



E a = c sP 2 -g s p{A a + A KQ ) 



Es = c s p\ - g s p 3 A K3 , 



E n 



CpPl- 



(61) 
(62) 

(63) 

(64) 

(65) 
(66) 



is given in Table V and will be discussed later. As far 
as the isospin dependence is concerned, the absence of 
non linearities (or couplings) in isovcctor fields leads to a 
rather poor isospin dependence since B 3 is constant. The 
introduction of a coupling of the isoscalar fields with the 
isovector p field corrects this fact. Finally, the main dif- 
ference is again in the functional describing the effective 
mass but now not only the expansions are made for dif- 
ferent quantities, the mass in relativistic approach and 
the inverse mass in the Skyrme model, but also the rel- 
ativistic approaches lead to a much richer functional of 
both pi and Tj. 



Then grouping the terms Ki, Eu and Eg leads to the 
following simplification 

Y J K l + E a +E s = (Mpi+Ti) +E' a + El 



t—p,n 



with 



T = 

K = 

E' s = 



2M* 



-c s p 2 + c S3 p a + c Si p A , 
-cspI, 



where the effective mass is given by equation (38) with 
(58) and (59). 

The non-linear effects are included in the additional 
parameters c S3 — Anc^/igl = ngl/Qm & s and c S4 = 
2Xci/3g 4 s = \gil2Aml 

The binding energy density B = E — Mp, reads now 



C. Non-linear ap and cop couplings 



Still a different model includes non-linear a — p and 
uj — p couplings [23, 24] which allow to change the density 
dependence of the symmetry energy. In the correspond- 
ing Lagrangian a new coupling term C aup is added: 



9% 



nKgtr + Kg 2 v v^]. (70) 



We have followed the prescription of [23] so that the cou- 
pling Aj is chosen in such a way that for kp — 1.15 fm _1 
(not the saturation point) the symmetry energy is 25.68 
MeV like in the NL3 paramctrization. We start by set- 
ting A s = as in [30]. In this case, the equations of 
motion are not the standard ones, once two of them be- 
come coupled and, for this reason, they are reproduced 
as follows: 



B 



non— rcl. 



EL + E„, + e' x + E n 



(67) 



and we get for B\ and B 3 



B 1 (p) = (c v -c s )p 2 + c S3 p 3 + (c S4 + c V4 )p 4 , (68) 



B 3 (P) 



-c$. 



(69) 



As in the simplest QHD model, taking advantage of the 
non-relativistic limit we have expressed the scalar density 
and the relativistic kinetic energy density as functionals 
of pi and Tj. Then the energy also becomes a functional 
of pi and Tj of the form of a Skyrme functional. The in- 
teraction part presents a very similar structure. In Figs. 
4a), b) and c) we show how the approximation works 
for NL3, TM1 and NL<5. While for TM1 the present ap- 
proximation is good, for NL3 and NL<5 we would have 
to include higher orders in the density expansion to im- 
prove the results. The non-linear coupling in the a and 
ui fields have introduced higher order terms in the poten- 
tial energy of symmetric matter. However because of the 
perturbative approach we have taken we are restricted 
to a polynomial density dependence. To go beyond this 
limitation, we have also fitted Skyrme parameters on the 
exact potential energy functional. The results of this fit 



9vV 
9p 



9l 



^L[p-2g v Vog 2 p b 2 A v ] 



bo = 



m 



Am 2 p 



[p 3 -4g p b glV 2 A v ] 



If the A s is not assumed to be zero but no non-linearity 
is taken into account either in the a or in the oj field (i.e. 
k = 0, A = and £ = 0), the three fields are coupled in 
the following way: 



g s <Po 

g v V 

9p, 

TT°0 



JhL[p-2g v V g 2 p b 2 A v 



[p s - 2g s (t>og 2 p b 2 A 



m 



^2 [(% - 4 -9pV5 s 2 0o A, + g 2 v V 2 A v )] 



This set of equations should be solved self consistently. 
However, if the A s and A„ are small we can solve the 
problem perturbatively introducing </> = -^%-p s , Vq = 



g v 



p and bo 



9p 



p 3 in the right hand side of the above 
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equations 



D. Density dependent coupling parameters 



g s (j>o = 2c s p s [1 - Mc s c 2 p plK s ] , 

g v V = 2c v p [l - 64c v c 2 p pjA v ] , 

^b = 2c pP3 [1 - Mc p {c 2 s pl A s + c 2 vP 2 A v )] . 

then the binding energy reads 



B 



E 

i— p,n 



Ti + E' a + E u + E' 5 +E P + E psv , (71) 



where the last term is the interaction energy between the 
p and the a and w fields 

Epsv = g 2 p b 2 [A s g 2 4>l + A v g 2 v V 2 }. 



In this case, 








n 

2Mf 


(72) 


K 


= -CsP 2 , 


(73) 


E u 


= c v p 2 - 2c vp p 2 pl, 


(74) 


E' s 


= -csPl, 


(75) 


E P 


= c pP\ - 2c v P P 2 pI - 2c sp p 2 pl, 


(76) 


E 


= c sp p 2 pj + c vp p 2 p 2 3l 


(77) 



with c vp — QAc 2 c 2 p A v , c sp = QAc 2 p c 2 A s . The effective mass 
can again be approximated by the leading terms of equa- 
tion (38). It should be noticed that the last term of E psv 
partly cancels the E p correction leading to 

E p + E psv = c p p\ - c vp p 2 p\ - c sp p 2 p\. 

The B coefficients are then written as 

= (c v -c s )p 2 + c S3 p 3 + c S4 p 4 (78) 
b 3(p) = Cp-cs- (3c vp + c sp )p 2 , (79) 

where we also include the contributions of the non-linear 
a terms. These expressions increase more rapidly with 
density than the corresponding Skyrme functional con- 
tribution given by equations (5-6): for the relativistic 
model a term with a a — 2 exponent would be neces- 
sary while the usual range of the Skyrme parameter a is 
below 1 in order to not present a too strong incompress- 
ibility. In Fig. 4 we see that the present approximation 
works well for B\ and the effective mass M* but fails to 
give a reasonable description of the binding energy. This 
was expected because NLojp is just NL3 with non linear 
Lup terms. The limitation of the present approximation is 
also clear for the B% term shown in Fig. 6a). Since we are 
looking at a low density expansion we have also directly 
fitted Skyrme parameters on the relativistic symmetry 
energy, Table V. 



Next we consider two models with density dependent 
coupling parameters, respectively the TW model [13] and 
the DDFM [10, 14] which also includes the 5 meson. These 
two models do not include self-interaction terms for the 
meson fields (i.e. k = 0, A = and £ = ). The only 
difference comes from the replacement of g coupling con- 
stants for the density dependent coupling parameters T s , 
T v , T p and Ts which are adjusted in order to reproduce 
some of the nuclear matter bulk properties, using the 
following parametrization for the TW model 



with 



and 



T t (p) =T t (p sat )h l (x), x = p/p sat , (80) 



, , , l + hix + di) 2 . 



h p (x) = cxp[-a p (.T - 1)], 



(82) 



and Ts(p) = 0, with the values of the parameters m^, T i: 
en, bi, Ci and d i7 i — s,v,p given in [13]. For the DDrM 
model we consider the TW parametrizations of T s and 
r„ and for the other two mesons we take 



hi(x) = en ex.-p[-b t (x - 1)] - c t (x - di), 



p, 8. 



Such density dependences in the coupling parameters do 
not affect the expression for the energy functional but 
of course affect its derivative such as the pressure or the 
chemical potentials. The latter ones are given by 



Pi = Vi + T V V + r l3 ^-b Q + £* 



(83) 



where the rearrangement term is 
dT v 



dp 



0T o b Q 8T S 8T S 

P3 7T 7T- Ps 00 TT~ Ps3 °3- 



dp 



dp 



dp 



As already discussed in the literature [10, 11, 15], the 
rearrangement term is crucial in obtaining different be- 
haviors in many quantities related to the chemical po- 
tentials or to their derivatives with respect to the den- 
sity as compared with the more common NL3 or TM1 
parametrizations . 

The binding energy functional is given by equation (67) 
with 



B^p) = (C v -C s )p 2 = C 0P 2 , (84) 

B 3 (p) - C p -Cs = D , (85) 

M* = M-2C s p-r 3i 2C SP3 , (86) 
= M-2C eS p-T 3i 2D cfiP3 , 

where the Ci coefficients are computed replacing the cou- 
pling constants gt by fj (p). 
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FIG. 3: Coefficients C; and A for (a) the TW model and (b) 
the DDH5 model. 



In Figs. 5 we show the binding energy, the coeffi- 
cient Bi(p) and the effective mass obtained in the non- 
relativistic approximation and the corresponding exact 
values. We see that the low density non-relativistic ap- 
proximation works quite well both with and without the 
5 meson. The same is true for the B3 contribution. In 
Fig. 6b) the exact and approximate coefficients B^(p) 
are plotted for these density dependent models (TW and 
DDH<5). 



E. Comparison between relativistic and Skyrme 
functionals 

We now compare the non-relativistic functional ob- 
tained from the relativistic models described so far with 
the Skyrme functional. 

In Table III we have collected the terms B\, B 3} Mi 
and M3 for the models we have considered in the present 
work. We notice that for the Skyrme forces the inverse of 
the effective mass is parametrized according to (7) while 
for the relativistic models we take a similar expression 
for the effective mass, Eq. (38). This fact explains the 
minus sign difference between the two types of models for 
the Mi and M3 columns. Except for the Walecka model, 
all models have for the isoscalar interaction contribution 
a parabolic term plus a higher order term on the den- 
sity. This second contribution appears explicitly through 
p 3 and p 4 terms, or implicitly through the density de- 
pendence of the coupling parameters for TW and DDFM. 
The isovector interaction contribution has a much poorer 
parametrization in the relativistic models: B 3 is gener- 
ally constant because D3 = 0, except for the NLojp model 
and again the TW and DDH5 models due to the density 



dependence of the coupling parameters. 

In Table IV we show the values of these coefficients for 
the models that we have discussed. One can see that Co 
and C e ff are of the same order of the corresponding pa- 
rameters of some Skyrme models shown in [3, 29]. The 
C3 coefficients are normally twice as large for relativistic 
than for non-relativistic models and the D coefficients 
for relativistic models are half of the coefficients of the 
non-relativistic ones. An immediate conclusion already 
referred is the poor parametrization of the isovector chan- 
nel in the relativistic models: for most models both D 3 
and D e ff are zero. Some comments with respect to the 
isoscalar channel are also in order: for the relativistic 
models the scalar kinetic contribution, denned by C e s is 
higher. This is due to the smaller effective mass within 
these models. The saturation is possible with an over- 
all larger binding for the isoscalar channel. This channel 
has an attractive term from the two-body force and a re- 
pulsive three-body (n-body) contribution. A larger bind- 
ing may be obtained with a stronger two-body attractive 
potential (Co) or a weaker three-body contribution (C3 
together with the a exponent). 

The coefficients of the density dependent models dis- 
cussed before could not be included in the above table 
since they are not fixed quantities. In Figs. 3 the coeffi- 
cients for these models, defined in equations (84-86), are 
plotted. Their values agree with the ones already given 
in Table IV for density values p/po > 0.2. 

In order to better compare the relativistic models, 
and in particular the parametrizations with density de- 
pendent coupling parameters, with several Skyrme force 
models we have also fitted the exact B\ and B% by the (5) 
and (6) expressions in the density range - 0.1 frn -3 us- 
ing, whenever appropriate, the same value for a in both 
expressions. In Table V we give the results of these fits. 

Some conclusions can be drawn from the values in Ta- 
bles IV and V. First of all we conclude that the isoscalar 
channel in relativistic mean field models has a quite com- 
plicated density dependence, both the interaction and 
momentum dependent term. In Table IV we only give 
the coefficients of the first terms of the expansion, which, 
as discussed before, works quite well at low subsatura- 
tion densities for QHD-II and TM1 but not so well for 
NL3, NL<5 and NLwp. Using expressions (5) and (6) to 
parametrize the isoscalar and isovector interaction term 
we have obtained for all the relativistic models a non- 
integer coefficient a smaller than 1 except for the QHD- 
II. This is a special case which, as we know, does not 
describe correctly nuclear matter properties, namely, it 
predicts a very large compressibility. All other a values 
are smaller than 1 but not so small as the correspond- 
ing parameter in Skyrme forces with good performances 
which are generally below 0.2. This small value of a in 
Skyrme forces controls the compressibility and is gener- 
ally taken equal to 1/6. Another term which has a very 
systematic behavior is the isoscalar momentum depen- 
dent term described by C c ff: this coefficient, except for 
the NL(5, is larger than 50 and maybe as large as ~ 100 
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MeV fm 5 for the TW. This occurs at low densities and is 
compensated by an extra binding coming from de attrac- 
tive term described by Cq. Non-rclativistic models have 
a similar behaviour, i.e., C e s ~ 50 MeV fm 5 . It should 
be noticed however that TW has a richer density de- 
pendence which is parametrized by the coefficient C e ff,3 
not present in the Skyrmc forces. On the other hand, 
relativistic models have generally a very simple isovector 
channel. The inclusion of the <5-meson brings in the extra 
degree of freedom missing in the momentum dependent 
terms but not in the interaction term if only linear terms 
are included for the mesons. Models with density depen- 
dent coefficients such as the TW include automatically 
a larger density dependence in the isovector term of the 
interaction. When we compare D3 for TW with D3 for 
SLy230a, NRAPR or LNS we observe a similar behavior. 

VI. CONCLUSIONS 

In the present work, we have compared relativistic 
mean-field nuclear models the relativistic models with 
the non-relativistic Skyrme models. We have shown that 
for the isoscalar channel the relativistic models behave 
in a similar way, and generally different from the non- 
relativistic description. This is true for the binding en- 
ergy, isoscalar interaction term and effective mass. The 
relativistic density dependent models give the closer de- 
scription to the one obtained by the non-relativistic mod- 
els. 

The isovector channel has proved to be a different prob- 
lem: there is a quite big discrepancy even between models 



within the same framework. This is the least known part 
of the nuclear interaction which we expect to determine 
with the future radioactive beams. Relativistic models 
have generally a very poor description of the this chan- 
nel. 

We have next tried to obtain a low density expansion 
of the relativistic models with a parametrization similar 
to the one used for the Skyrme interactions. The energy 
functional of the relativistic models depends not only on 
the isoscalar and isovector particle and kinetic densities 
but also on the isoscalar and isovector scalar densities. 
Only in the low density regime these densities reduce to 
the respective particle densites. We have shown that for 
some models already in the subsaturation density expan- 
sions it is necessary to include many terms in order to 
get good agreement. For the low density range for which 
there is a good agreement between the low density expan- 
sion and the exact values, we have shown that some of 
the coefficients of the relativistic and the non-relativistic 
models are of the same order of magnitude. However 
we have shown that for the first ones the scalar kinetic 
contribution is higher. The saturation is possible with 
an overall larger binding for the isoscalar channel which 
is due to a stronger two-body attractive potential or a 
weaker three-body contribution. 
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FIG. 4: Comparison between several exact and approximate 
physical quantities: a) binding energy density, b) -Bi(p) coef- 
ficient and c) effective mass M* for relativistic models with 
constant couplings. 



Model B/A po K £ sym M*/M 

(MeV) (fm~ 3 ) (MeV) (MeV) 

SIII [3] 15.851 0.145 355.5 28.16 0.76 

SLy230a [3] 15.989 0.16 229.87 31.97 0.697 

NRAPR [4] 15.86 0.16 225.7 32.79 0.7 

LNS [25] 15.32 0.175 210.85 33.4 0.825 

TABLE I: Nuclear matter properties of the Skyrme forces 
used in the present work 



Model 


B/A 


pa 


K 




M*/M 




(MeV) (fm" 3 ) 


(MeV) (MeV) 




Walecka [19] 


15.75 


0.192 


540 


22.1 


0.556 


NL3 [20] 


16.3 


0.148 


272 


37.4 


0.60 


TM1 [21] 


16.3 


0.145 


281 


36.9 


0.63 


NL5 [22] 


16.0 


0.160 


240 


30.5 


0.60 


TW [13] 


16.3 


0.153 


240 


32.0 


0.56 


DDH<5 [14] 


16.3 


0.153 


240 


25.1 


0.56 



TABLE II: Nuclear matter properties of the relativistic mod- 
els used in the present work 



TABLE III: Expressions for the coefficients Bi(p), B- 3 (p), 
Mi(p) and M-z(p) . For the Skyrme parametrization we take 
(7) which refers to M* _1 . 



model 




Bi(p) 




B 3 (p) 


M 1 (p)/p 


M 3 (p) 


Skyrme 


C p 2 


+ c 3P °+ 2 




Do + D z p a 


2C c ff 


2D cS 


QHD-II 


(c„ 


~Cs)p 2 




c p 


-2c s 




NL3, TM1 (c„ 


- c s )p 2 + c 


s3P 3 + {Cs4 


+ c vi )p 4 


Cp 


-2c s 




NL5 (c„ 


- Cs)p 2 + C 


s3P 3 + (c s 4 


+ c vi )p 4 


Cp — Cg 


-2c 3 


-2c s 


NLojp, NLo-cj 


(c v — c 3 )p 


+ c s3 p 3 + 


CsAP Cp 


(3cup + c S p)p 


-2c s 




TW 


(C v 


- C s )p 2 




c P 


-2C 3 




DDH<5 


(C v 


- Cs) P 2 




C p — Cs 


-2C S 


-2C S 
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TABLE IV: Coefficients d and Di obtained in the present 
work and results from several Skyrme models. 



mnrlpl 


^0 


c 3 


D 


Dq 


(7 «■ 

^eff 


D a 
^cfr 


(j 




CMeV fm 3 ) 


(MeV fm 3+3<T ) (MeV fm 3 ) 


(MeV fm 3+3 °") 


(MeV fm 5 ) 


(MeV fm 5 ) 




OHr> TT 


ouo.oo 


n 
u 


oy.ou 


n 


Q7 


n 

u 


\ 


NT ^ 


-O 1 1 .04: 


2482.69 


131.44 


n 

u 


R8 01 

UO.Ul 


n 

u 


i 
i 


NT,m/i A — 01 


-511.34 


2482.69 


149.22 


-5975 74 (n — 2) 


68.01 





1 


NLwp,A„ = 0.025 


-511.34 


2482.69 


189.76 


-24161.25 (cr = 2) 


68.01 





1 


TM1 


-479.21 


1571.51 


138.34 





64.90 





1 


NL<5 


-482.07 


2369.38 


63.82 





44.64 


10.81 


1 


SIII [3] 


-426.28 


875 


268.08 





44.38 


-30.63 


1 


Ski' [29] 


-396.49 


903.97 


208.42 





12.98 


-20.99 


1 


SLy230a [3] 


-933.84 


862.69 


1015.89 


-1392.89 


56.37 


56.37 


1/6 


NRAPR [4] 


-1019.89 


940.13 


449.80 


-398.68 


57.015 


-27.992 


0.14416 


LNS [25] 


-931.86 


911.76 


349.62 


-283.17 


25.05 


-19.5 


1/6 



TABLE V: Coefficients obtained from the fitting to the exact 
Bi and B 3 expressions 



model 


Co 


c 3 


Do 


D 3 


C e ff,0 


C c ff,3 


D cB 


a 




(MeV fm 3 ) (MeV fm 3+3<T ) 


(MeV fm 3 ) (MeV fm 3+3<T ) 


(MeV fm 5 ) (MeV fm 3+3<T ) (MeV fm 5 ) 




QHD-II 


-297.0 


253.625 


59.30 





50.97 








2.5 


NL3 


-625.25 


546.28 


131.44 





68.01 








0.318 


TM1 


-482.84 


546.07 


138.34 





64.90 








0.606 


NL<5 


-594.83 


802.09 


63.82 





44.64 





10.81 


0.536 


TW 


-486.09 


613.50 


251.92 


-746.72 


96.29 


-136.80 





0.767 




FIG. 5: Comparison between several exact and approximate 
physical quantities: a) binding energy density, b) Bi (p) coef- 
ficient and c) effective mass M* for relativistic models with 
density dependent couplings. 




FIG. 6: Bz(p) coefficient for a) a non-linear coupling model 
(NLwp) and b) two density dependent models TW and DDH<5. 



